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ABSTRACT 

The reionization of the intergalactic medium (IGM) was likely inhomogeneous and extended. 
By heating the IGM and photo-evaporating gas from the outskirts of galaxies, this process can 
have a dramatic impact on the growth of structures. Using a suite of spherically-symmetric 
collapse simulations spanning a large parameter space, we study the impact of an ionizing 
ultraviolet background (UVB) on the condensation of baryons onto dark matter halos. We 
present an expression for the halo baryon fraction, /b, which is an explicit function of: (/) 
halo mass; (ii) UVB intensity; (z'zz) redshift; (z'v) redshift at which the halo was exposed to 
a UVB. We also present a corresponding expression for the characteristic or critical mass, 
M cr it, defined as the halo mass which retains half of its baryons compared to the global 
value. Since our results are general and physically-motivated, they can be broadly applied to 
inhomogeneous reionization models. 



Key words: 

evolution 



cosmology: theory - early Universe - galaxies: formation - high-redshift 



> 

t> 

o 

m 

> 

X 
C3 



1 INTRODUCTION 

As the first galaxies formed, their ionizing UV radiation carved- 
out HII regions from the neutral intergalactic medium (IGM). As 
structures continued to form, these HII regions grew and over- 
lapped, eventually permeating the entire Universe. This global 
phase change, known as reionization, is expected to be fairly ex- 
tended and inhomogeneous (e.g. Furl anetto et aL 2004; Zah n et aT| 
|201 1[ >. The resulting ionizing ultraviolet background (UVB) heated 
the IGM to temperatures of ~ 10 4 K, photo-evaporating gas out 
of shallow potential wells and affecting its cooling properties (e.g. 



Shapir oeFai.|1994| |Miralda-Escude & Rees|1994{ |Hui & Gnedm] 



1997| >. Therefore the gas reservoir available for star formation is de- 
creased for low-mass galaxies in the ionized IGM (e.g. |Couchman| 
|&Rees|1986l|Rees|1986t|Efstathiou|1992) . 

This UVB feedback mechanism can have several important ef- 
fects: ((') by suppressing star-formation inside cosmic HII regions, it 
results in a more uniform reionization morphology (e.g. McQuinn 
|et al.|200 7 ); (ii) it can delay the completion of reionization, helping 
to alleviate potential tension between a reionization with a mid- 
point at^~^10j&mSs^^L]2^n i, and an end-point at z ~6-7 
(e.g. |Mortlock et al. 2011 Pentericci et al.||201 1( |Dijkstra et al.| 



201ll|Bolton et al.|2011||Schroeder et al.|2013||Ono et al.|2012 



explaining the apparent dearth of observed satellite galaxies in the 
Milky Way (e.g. |Klypin et al.|[l999l |Moore et al.1[T999"l |Bullock| 
|et al.| [2000 Somerville 2002J and dwarf galaxies in the field (as 
inferred from the HI ALFALFA survey; [Papastergis et al |201 1| 
|Ferrero et al.|2012) . 

Since semi-analytic calculations (e.g. |Shapiro et al.| |1994| 
|Gnedin & Hui|19 98 ) were used to motivate a suppression of star 
formation in galaxies with M < 1O 9 M0, several attempts have 
been made to further quantify this effect, both using spherically 
symmetric simulations iThoul & Weinberg 1995; Kitayam a et al.| 



(Hi) it can suppress the baryon content of local dwarf galaxies, 



|2000||Dijkstra et al.|2004) and three-dimensional cosmological hy- 
drodynamic simulations (e.g. |Quinn et al.||1996| | Weinberg et al.| 
\TWf\ [Navarro & Steinmetz|1997||Gnedin|2000||Hoeft et al.|2006T 
IQkamoto et al.|2008| >. |Gne din (2000) (see also |Naoz et al.|2009| > 
used numerical simulations to argue that the resulting decrease of 
baryonic content of halos was well predicted by the linear theory 
filtering mass (Gnedin & Hui 1998). However, subsequent simu- 
lations showed that this prescription should to be revised in the 
low-redshift, post-reionization regime (Hoeft et al. 2006] |Okam6tol 
|et al.|2 008 1, taking into account the cooling properties of gas at the 
outskirts of galaxies. 

Unfortunately, quantifying UVB feedback is challenging. Semi- 
analytic approaches do not include all of the relevant physics (e.g. 
non-linear growth of perturbations, cooling processes of the gas). 



* email: emanuele.sobacchi@sns.it 

1 Note that extended reionization scenarios could also result from feed- 
back from early sources of X-rays I Ricotti & Ostriker 2004 , Mesinger et al. 



2012), or from reionization "stalling" when the HII regions grow to sur- 
pass the attenuation length of ionizing UV photons ( Furlanetto & Mesinger 
[200^|[Croc7ani et al.|2011||Alvarez & Abel|2012| . 
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M = 1.2xl0*M o atz = 6 




Figure 1. Evolution of the gas shell enclosing half the baryonic mass of 
a 1.2 X 10 9 Mq halo collapsing at z = 6, for several of our runs. The 
radius is expressed in units of i?M AX : the maximum radius reached by the 
outermost shell (enclosing the total halo mass at z = 6) in the run without 
the UVB. Similarly, the time is expressed in units of t co ii> the cosmic time 
at the collapse redshift (z = 6) of the outermost shell. Dotted: no UVB. 
Dot-dashed: the UVB (with J21 = 1) is turned on at jzjn = 9. Short- 
dashed: the UVB (with J21 = 0.01) is turned on at zm = 11. Dashed: 
the UVB (with J21 = 1) is turned on at zm =11. Long-dashed: the UVB 
(with J21 = 1) is turned on at Zm = 18. 

On the other hand, cosmological simulations often focus on a 
particular reionization history, fixing the UVB evolution and/or 
the reionization redshift. Our poor understanding of the details of 
reionization motivates a flexible approach that is not dependent on 
a particular model, and can therefore be broadly applied. 

As an intermediate approach, in this paper we run suites of fast, 
ID cosmological collapse simulations, exploring a wide parameter 
space motivated by the inhomogeneity of reionization. Fitting these 
results, we present an expression for the halo baryon fraction, /b, 
as well as the critical mass M cr it, defined as the halo mass which 
retains half of its baryons compared to the global value. 

In this work, we focus on halos with masses of 10 8 < M/Mq < 
10 10 , which are relevant in the advanced stages of reionization (e.g. 
|Haiman et al.|1997) . These halos are massive enough to host gas 
collapsing via the atomic cooling channel, and yet small enough to 
be susceptible to UVB feedback. Furthermore, we do not attempt to 
model baryonic suppression through internal feedback processes, 
such as internal radiation or winds (e.g. Springel & Hernquis t|2003| 
Hopkins et al. 201 1 ). These additional processes require parameter 
studies of detailed interstellar medium physics of high-2 galaxies; 
isolating their contribution is best done through a separate line of 
investigation. 

The paper is organized as follows. In S|2]we describe the collapse 
simulations we use. In *j5]we use results from these simulations to 
fit formulae for M cr it and /b, also highlighting the limitations of 
our approach. Finally, in S|4]we present the conclusions of the work. 
Throughout we assume a flat ACDM cosmology with parameters 
(O m , Ov, £^b, h, cr 8 , n) = ( 0.27, 0.73, 0.046, 0.7 0.82, 0.96), con- 
sistent with WMAP results (Kom atsu et al.|201l) . 



2 HYDRODYNAMICAL COLLAPSE SIMULATIONS 

We use a ID simulation to study collapsing, cosmological perturba- 
tions. The code evolves a mixture of dark matter (DM) and baryon 
fluids by moving concentric spherical shells of fixed mass in the 



radial direction. It includes cooling through excitation, ionization, 
recombination and free-free emission and by Compton scattering 
with CMB photons (relevant at z > 7; |Dijkstra et al.|2004| . For 
more details on the code please see |Thoul & Weinberg|<|1995|>. 
As in previous works (e.g. [Thoul & Weinberg||1995| |Dijki 



g| |iw5} [Dijkstra 

|et al.|2 004 , Mesinger & Dijkstra 2008} we start wi th typical pro- 
file around a 2-a pealj^jin a Gaussian random field ( jBardeen et al 



[1986) , sampled by 6000 (1000) shells for the DM (baryons). To set- 
up the runs, we fix the total halo mass M collapsing at redshift z. 
Note that due to the hierarchical nature of structure formation, the 
halo mass increases with time, so M at 2 corresponds to a smaller 
Mi < M at zi > z. [J] We include an isotropic UVB instanta- 
neously turned on at redshift zm and parametrized as 



J(v) = J21 {v/va)~ a x 10" 



erg s 



(1) 



McGreer 



where is the Lyman limit frequency and a = 5 corresponds to 
a stellar-driven UV spectrum (e.g. |Thoul & Wei nberg 1996}. We 
consider a range of 2in= 9-16, as well as three different intensi- 
ties: J21 = 0.01, J21 =0.1 and J21 = 1, spanning the expected 
range of interest]^] This parameter space is motivated by paradigm 
of inhomogeneous reionization: as reionization proceeds, ioniza- 
tion fronts pass through different regions of the IGM at different 
redshifts, zm, until they finally overlap at z >5-6 (e.g. 
|etal.|2011| >. 

The effect of the ionizing background on the evolution of a gas 
shell is shown in Fig. [T] This shell encloses half of the baryonic 
mass of a M = 1.2 x 10 9 M Q halo collapsing at 2 — 6; this total 
halo mass corresponds to, e.g., M = 1.2 x 10 s Mq at t = 3£ C oll/4 
(2 ~ 7.4) and M = 2xlO 7 M at* = t coU /2(z ~ 10). The dotted 
curve shows the evolution of the shell with no UVB. As expected 
in the (pressureless) spherical collapse scenario (e.g. |Gunn & Gott] 
|1972| l, the maximum radius of this shell is reached at half of its 
collapse time. 

The dot-dashed curve shows the evolution of the shell when a 
J21 = 1 UVB is turned on at zm = 9. The gas is heated by the 
UVB and the resulting pressure support delays the collapse. Inter- 
estingly this shell (enclosing half of the baryonic mass) collapses 
at t ~ tcoli, indicating that M = 1.2 x 1O 9 M0 corresponds to the 
so-called characteristic or critical mass (see below) for this param- 
eter combination: 2, zm and J21. Similar behavior is shown with 
the short-dashed curve, corresponding to the run when the back- 
ground is turned on earlier, zm = 11, but with a reduced intensity, 
J21 = 0.01: the evolution is delayed since earlier times but the 
pressure support is lower, therefore this shell also collapses at t ~ 



2 We find that our results are not very sensitive to the choice of the peak 
height. In particular, switching to a 1.5 (2.5) sigma peak lowers (raises) by 
~ 15% the derived value of M CT n. 

3 To facilitate comparisons below, we define the halo mass M indepen- 
dently of the UVB. Since the accretion of DM (making up most of the 
mass) on halo scales is not sensitive to the UVB, this convention effectively 
translates to M = (S7m/fiDM).MDM- 

4 Measurements of the Lya forest imply J21 ~ 1-0.1 at z ~2-6 (e.g. 
|Bolton & Haehnelt|2007"l|Calverley et al.|201 1| [Haardt & Madau|2012fr . 
One might expect this decreasing trend of the UVB to continue to higher 
redshifts, since there are fewer collapsed structures and the mean free path 
of ionizing photons is decreasing with redshift (e.g. McQuinn et al. 201 1 
|Cowie et al.|2 01 1 1. However the drop in the number of galaxies could be 
compensated by an increase in their ionizing efficiency, due to e.g. a higher 
escape fraction or a top-heavy stellar initial mass function (e.g. | Schaerer| 
|2002||Wise & Cen|2009[|Yajima et al.|201 l|[Ferrara & Loeb|2012) . Since 
the evolution of the UVB intensity is very uncertain at z > 6, we explore a 
large range of intensities. 
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t C oii. These two points in our parameter space, (M, z, Zm, J21) = 
(1.2 x 1O 9 M , 6, 9, 1) and (1.2 x 1O 9 M , 6, 11, 0.01), both result 
in a baryon content which is ~ 1/2 of the value expected without 
aUVB, i.e. f b ~ 1/2. 

The medium-dashed curve shows the evolution of the shell when 
J21 = 1 and zm = 11: the collapse is more significantly delayed 
and the corresponding halo baryon fraction is / D < 1/2. This is 
also true for the long-dashed curve, corresponding to the same in- 
tensity J21 = 1 but a higher zm = 18. In this case the shell ex- 
pands slightly after exposure to the UVB and is pressure-supported, 
with the radius evolving very slowly. However, as the enclosed po- 
tential well deepens from further DM accretion, the shell eventually 
starts to collapse again at t ~ t co ii- 

The evolution is simpler in smaller halos with M ~ 1O 8 M0 at 
z = 6 (not shown in the figure) since the gas can be more easily 
photo-evaporated out of the shallower potential wells. Shells en- 
closing even half of the mass "peal-off" at z = Zm, effectively 
escaping back into the IGM. 



3 RESULTS 

Armed with our suite of collapse runs, we now present the 
two main results of this paper: a formula for the critical mass 
M cr it {J>z, Zm) (defined as the halo mass retaining half of its 
baryons; §3.1} , and its generalization to the halo baryon fraction 
fb (M, J, z, zm) (S 3.2 1. Our approach allows us to explicitly in- 



Mindful of these trends, we assume the functional form: 



elude dependences of these quantities on J and zm- This is im- 
portant both for its generality, and also since all reasonable reion- 
ization models are very inhomogeneous (e.g. Mesinger et al. 2012 
|Sobacchi & Mesinger|20T3) . 

3.1 Critical Mass 

As an ionization front sweeps through a patch of the IGM at 
z = zm, the temperature of the gas is expected to jump from 



lz>z IN 



0, to T. 



z<z IN 



10* K (e.g. 



Barkana & Loeb 2001 



roughly within a sound-crossing time.Using linear theory in such 
a scenario, |Gnedin & Hui ( 19981 (Appendix A) motivate the sup- 
pression of gas on scales smaller than the "filtering" wavenumber: 



3 

Wkj 



1 + 4 



l + z 
1 + zm 



5/2 



- 5 



l + z 
1 + zm 



(2) 



Gnedin 



2000 



where kj corresponds to the classical (instantaneous) Jeans 
wavenumber at T Z<ZIN . Although there is no obvious motivation 
why this linear theory derivation should extend to collapsing and 
cooling structures, the filtering scale was shown to agree well with 
hydrodynamical simulations at high redshifts (z > 6 
IQkamoto et al.|2008[|Nabz et al.|2009l >. 

We note from eq. |2](, that the temperature dependence (encoded 
in the Jeans wavenumber) factors out. Since the intensity of the 
UVB mainly affects the gas temperature, we consider a factorable 
expression M CI1 t(J, z, zm) = f(J)g(z, zm)- This is confirmed by 
our simulations a nd also consistent with former ones (e.g. Mesinger 



& Dijkstra 



2008 1, which motivate a power law form / (J) oc J 2i a . 
Furthermore, we expect that the asymptotic limit of M CT it(z <C 
zm) is approached independently of the precise value of zm- 
Therefore we can write Merit (J,z, zm) = J2i a gi ( 2 ) 92 (z, zm), 
with #2(2, zm) ~ 1 if 2 "C zm- Since the Jeans mass before 
reionization is much lower than the mass scales we are studying, 
we assume the additional limiting behaviour of gi(z, zm) — > at 
z — > zm- 



M cri t = M J 2 i a 



l + z 
10 



l + z 

1 + ZIN 



(3) 



where we treat Mo, a, b, c and d as fitting parameters ^ A \ 2 min- 
imization fitting to our simulation outputs results in: 

(M ,a,b,c,d) = (2.8 x 10 9 M Q , 0.17, -2.1, 2.0, 2.5) (4) 

Values of M cr it from our simulation runs are shown as points in 
Fig. [2] together with the curves corresponding the analytic expres- 
sion above (eq.[3}j4]l. The relative errors of this fitting formula with 
respect to the simulation values are all < 10%. We have tested the 
robustness of the fit by including additional data points with dif- 
ferent choices of z and zm (each spanning the range 7-19), and 
find that the best-fit parameters changed only within few percent 
of those in j4}. Note that our results are not very sensitive to the 
intensity of the UVB: increasing J21 by two orders of magnitude 
increases Merit only by a factor ~ 2-3. This weak dependence is 
expected given that the post-ionization collapse in this regime pro- 
ceeds isofhermally, and T z<ZlN ~ 10 4 K is only weakly dependent 
on J jDijkstra et al.|2004| >. 

The functional form for M cr it assumes a rapid transition at zm, 
i.e. r z>ZlN <C T z<ZlN . However, earlier populations of X-ray 
sources could have pre-heated the IGM, making the change in the 
Jeans mass at zm less dramatic. Of greater interest is the ability of 
the gas to cool efficiently, allowing it to continue contracting past 
the adiabatic stage, eventually forming stars. It is thought that the 
first galaxies formed in halos with virial temperatures T V1I <C 10 4 
K, governed by the molecular cooling channel. However, these ob- 
jects were very sensitive to highly uncertain feedback processes 
(e.g. |Machacek et al.||2003| |Kuhlen & Madau| [2005), making it 
difficult to estimate the effective minimum halo mass required to 
host galaxies, M m i n , at epochs around and before zm ■ Star-forming 
galaxies need to be massive enough to both: (i) host gas in the pres- 
ence of a UVB (M > Merit); and (ii) allow this gas to efficiently 
cool, collapse and form stars (M > M CO ol)- Therefore, a reason- 
able, general choice for M m i n would be: 



M„ 



max [M c , 



cool i -^^crit J ; 



(5) 



where M CO oi corresponds to the effective cooling threshold (includ- 
ing feedback effects). 

Our expression for M cr it is a good match to the "filtering" mass 
scale, as well as to numerical simulation at high-redshifts (z > 6; 
e.g. |Okamoto et al.||2008) . However our expression has a steeper 
low-redshift evolution with respect to the simulations of Okamoto 
|et al.| ( [2008] l, over-predicting the value at z = by a factor of w 10 
(see their appendix B). The particular realization of Okamoto et al. 
( 2008 1 prefers an M cr it which asymptotes to a constant virial tem- 
perature, i.e. Merit oc (1 + z) -3 / 2 . A possible reason for this dis- 
crepancy is the fact that our simulations are missing 3D substruc- 
tures, whose resulting potential wells would aid in the retention 



5 The same form with parameters (Mo, a, b, c, d) = (1.9 X 10 9 M Q , 0, 
—3/2, 3/2, 3) is obtained under the assumption that the baryons get pho to- 
evaporated over a sound crossing time at 10 4 K (c.f. |Shapiro et al. 2004 1 



200Myr 



M 



1/3 



l + z 
10 



0.15 



-1/3 



, 10 8 Mq , 

and that the halo grows over a longer time-scale. In this scenario M cr it is 
obtained solving the equation t ac (Af cr it) = *H ( z ) — *H (^in). where tjj 
is the cosmic time. 
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Figure 2. Critical mass versus redshift with different choices of 2tn; points show the output of the simulations and curves show the best fit of eq.[3]with three 
different values of J21: Dot-dashed: J21 = 0.01. Dashed: J21 = 0.1. Dotted: J21 = 1. 



of gas. Hence, we caution against over-interpreting eq. {3} at low- 
redshifts (z < 5), where halos increasingly grow by mergers and 
our simple, spherically-symmetric picture breaks down. In a com- 
panion paper \Sobacchi & Mesinger 2 013\ , we provide a general, 
analytic expression for the global average A-/ m i n (,z), considering a 
"fiducial " inhomogeneous reionization model and a constant T v - lT 
low-redshift limit. 



3.2 Halo Baryon Fraction 

We now generalize the expression for the critical mass, obtain- 
ing the halo baryon fraction, fb, normalized to the global mean, 
fib/^m. By definition, fb(M = A/ cr it) = 1/2, with the asymp- 
totic behaviour: f b (M < A/ cr it) = and f b (M > A/ crit ) = 1. 
We find our results are well-fit (with relative errors < 5%), by the 
following simple expression]^] 

/ b (M) = 2- AWM . (6) 

We note that this expression is a relatively steep function of M, thus 
justifying the common, step-function simplification of fb = — > 1 

at M cri t. 

In Fig. [3] we show the baryon fraction of halos with different 
masses and reionization histories. Points show the output of the 
simulation and curves show our fitting formula. We present a "fidu- 
cial" model, (z,zin, J21) = (7, 10,0.1), varying one parameter 
in each panel. In the left panel, we show the variation of fb with 
J21. As we have already discussed this dependence is weak. In the 
middle panel, we show the variation of fb with ztn- This depen- 
dence is stronger. In particular, regions only recently exposed to a 
UVB do not exhibit significant baryonic suppression. This is qual- 
itatively supportive of earlier claims (Mesinger & Dijkstra 2008) 
that UVB feedback does not significantly impact the progress of 
the advanced stages of reionization (note that this panel shows a 
"maximal" zm — z case of halos collapsing towards the very end 
of reionization, z = 7). These results also suggest that a potentially 
large scatter in the reionization redshift of halos (e.g. Alvarez et al. 
|2009| > may have an important impact on the evolution of their gas. 
Finally, in the right panel we show the variation of fb with z. The 
formation redshift has a very strong impact on the baryon fraction, 
consistent with claims that UVB feedback can be important in sup- 
pressing star formation post-reionization (e.g. Thoul & Weinberg 
[19961 |Benson et al.|2002| >. 



6 Our results are also well-fit with commonly-used expression suggested 
by Gnedin 1 2000 1 (his eq. 7, with a = 1). However, as that functional form 
is also not physically motivated, we use the simpler expression in eq. |6j. 



3.3 Limitations of our Method 

Using spherically-symmetric collapse simulations allows us to 
rapidly explore a broad parameter space, explicitly quantifying the 
dependence of the halo's baryon content on its reionization history. 
This is a critical improvement to previous work since reionization 
is likely to be very inhomogeneous. Our general approach facili- 
tates its inclusion in reionization models. Nevertheless, it suffers 
from several limitations: 

(i) Evolving background: Our parameter space only includes 
non-evolving UVB intensities. As mentioned above, one might ex- 
pect the ionizing background to decrease towards higher redshifts. 
Observations of the Lytv forest post-reionization imply a fairly con- 
stant UVB over the range 2 < z < 5 (e.g. Bolton & Haehnelt 
|2007| >. Estimating its evolution towards higher redsmlts is limited 
by our understanding of absorption systems and galaxy emissivity 
(e.g. |McQuinn et al.|2011| >. In any case, we show that feedback is 
fairly insensitive to the UVB intensity, so we do not expect any 
realistic evolution to affect our results. 

(ii) Self-shielding: Our models do not include self-shielding, 
which could decrease the amount of ionizing radiation reaching 
the inner parts of the structures, thus reducing the impact of UVB 
feedback (e.g. |Susa & Um emura 2004). Spherically symmetric 
models likely overestimate this effect since ionizing photons can 
likely penetrate inside structures following low HI column density 
paths. Furthermore, the UVB can strip away the outer layers of 
gas, rapidly exposing deeper layers which were previously self- 
shielded, with photo-evaporation occurring in roughly the same 
time-scale as expected in the optically-thin limit flliev et al.|2005| >. 
Regardless, we find that the evolution of gas at the outskirts of halos 
plays a dominant role in governing the UVB feedback, consistent 
with previous claims (e.g. Oka moto et al.|2 008 1. 

(iii) 3D structure: Although spherical symmetry is a decent 
approximation for early-collapsing structures (e.g. Bardeen et al. 
1986), this approximation becomes increasingly suspect towards 
lower redshifts. Substructures form earlier than the collapse of the 
analogous spherical shell, making it easier for them to retain gas. 
As discussed above, we suspect this to be an important consider- 
ation at 2 < 5 1 Okamoto et al. 2008). Further progress is likely 
to involve suites ot iU simulations (including radiative transfer), 
guided by the functional forms motivated above. 



4 CONCLUSIONS 

Using a spherically symmetric simulation, we study the baryonic 
content of atomically-cooled galaxies exposed to a UVB. Since the 
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Figure 3. The baryon fraction, /b, versus M for halos with different reionization histories. Points show the output of the simulations and curves show our 
fitting formula in eq. |6j. We choose a "fiducial" parameter combination: (z, zjn, J21) = (7, 10, 0.1) and show variation with J21, zin, Z in the left, middle 
and right panels (respectively). 



details of reionization are unclear, we explore a large parameter 
space of (i) halo mass; (ii) UVB intensity; (iii) redshift; (iv) red- 
shift of UVB exposure. The later is especially important given that 
reionization is likely very inhomogeneous. We present an analytic 
expression for the characteristic or critical mass M cr it, defined as 
the mass scale of halos retaining half their gas mass compared to 
the global mean. We also generalize these results, obtaining a sim- 
ple formula for the baryonic content of galaxies as a function of 
(i)-(iv) above. Our results can be readily applied to models of in- 
homogeneous reionization and high- z galaxy formation. 

We thank Mark Dijkstra for providing and technical assistance 
with the code we have used as well as for interesting discussions. 
We thank Matthew McQuinn, Stuart Wyithe, Zoltan Haiman, and 
Andrea Ferrara for the helpful comments on early drafts of this 
paper. 
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